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1 Introduction 



During the past two decades one of the challenging problem of QCD has been to understand 
theoretically and to observe experimentally a new nonperturbative regime - high parton density 
QCD. It has been argued in perturbative QCD (pQCD) that at low x ( high energies ) 
the density of parton increases in deep inelastic scattering ( DIS ) and reaches so high value 
that partons become densely populated in a hadron. Certainly, high density system of partons 
cannot be treated perturbatively but pQCD leads to a new scale ( mean transverse momentum of 
partons ) for high parton density QCD (hdQCD ) as well as to a hypothesis of parton saturation 
Intensive theoretical studies led to deeper understanding of physics of this system as well 
as to development of both perturbative and nonperturbative methods for such a system. 
However, only recently the first indications on parton saturation have appeared in HERA data on 
- behaviour of the F2-slope and on energy behaviour of inclusive diffractive dissociation in DIS 
0]. These data as well as their theoretical interpretation mark a new stage of our approach 
which needs more quantitive methods, than it was before, to make a reliable prediction for the 
experimental observations. 

This quantitive approach includes at least two steps: 



1. a derivation of the equation which will be valid in the full kinematic region; 



2. finding a solution to this equation. 



There are two versions of the nonlinear evolution equation at x — ^ 0. The first one was suggested 
as an obvious generalization of pQCD approach ||^ while the second has been obtained from 
semiclassical gluon field approach 0. Both of them describe correctly the DGLAP § evolution 
in the limit of low parton density as well as the GLR nonlinear evolution equation [|l| in the 
region of intermediate parton densities. They give the same limit of parton density saturation 
at very low values of x but they are different in particular description how system reaches this 
saturation. Yu. Kovchegov has recently proved ^ that the GLR equation itself (in the form of 
Eq.(2.18) in Ref. ) is able to describe the QCD evolution in the region of high parton density 
or, in other words, at very low values of x. His arguments, which we review briefly in the next 
section, are based on Mueller's idea that colour dipoles rather than quarks and gluons, are 
correct degrees of freedom at high energies ( low x ) in QCD. It was also shown that the solution 
to the equation, suggested in Ref. [0, coincides with the solution to the GLR equation. 

The goal of this paper is to solve the GLR equation, written in the form suggested in Ref. 0, 
in the full kinematic region including x — 0. The paper is organized as follows: In the next 
section we briefly discuss the nonlinear evolution equation for the high parton density QCD and 
we formulate our approach for searching a solution to this equation. In section 3 we introduce 
the new scale Qo{x) , which appears in hdQCD as a solution to the nonlinear equation, and we 
find the general solution to the equation for parton with > Qq{x), based on the approach 



developed in Ref. . Section 4 is devoted to a solution of the hdQCD equations in the region of 



very low x. We show that the parton density reaches the saturation limit and we find the explicit 
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analytic expression describing how the system approaches the saturation hmit. Our summary 
and conclusions are given in section 5. 



2 Nonlinear equation for high density parton system 
2.1 Kinematics, notations and definitions. 

In this subsection we introduce all notations and definitions as well as some kinematic relations 
that we will need in what follows. In general, we try to use the same notation as in Ref.p|. We 
would like to clarify the relations between different notations and definitions that have been used 
in this area of activity. This task is rather complicated since different notations have been used 
for the same physics observables, mostly due to different theoretical background of authors. In 
this paper we explore the main physical idea of Ref . |]T0| , namely, that correct degrees of freedom 
in QCD at low x ( high energy ) are colour dipoles rather than quarks and gluons that explicitly 
written in the QCD Lagrangian. This statement means that a QCD interactions at high energy 
do not change the size and the energy of a colour dipole. Therefore, the majority of our variables 
and observables is related to distributions and interactions of the colour dipoles in a hadron. 



1. Xik = |xi — Xk| is the size of the dipole which consists of a quark "f at Xi and an antiquark 
"fc" at Xk; OT, in other words, it is the transverse separation between quark and antiquark 
in a colour dipole; 

2. r = In ^ = In^, where ( Ql) is the transverse momentum of quark in the dipole 
with size x ( xq ) respectively; 

3. X is the Bjorken variable for a dipole , x = Q'^/W'^, where W is the dipole energy; 

4. y = ln(a;o/a;), where xq is defined for us the region of low x. We consider that for x < xo 
all the typical features of low x physics should be seen. Practically, xq ~ 10^^; 

5. xG{x, Q"^) is the gluon structure function for a nucleon; 

6. xGa{x, Q^) is the gluon structure function for a nucleus with A nucleons; 

7. bt is the impact parameter for the dipole scattering; 

8. The main physical observable, which we are going to discuss in this paper, is the density of 
dipoles ( A^(x, bt, y)) with the size x and energy x [y ) aX the impact parameter bt, 
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9. In the parton approach ( parton = dipole ) the dipole density N{-x,bt,y) is related to 
the dipole scattering amplitude if we assume that this amplitude is dominantly imaginary. 
Therefore, A^(x, bt,y) can be normalized using the unitarity constrains: 

2N{^,h,y) = N\^,h,y) + G*"(x, 6^, 2/) , (2.1) 

where G*" is the contribution of the inelastic processes to the scattering of a dipole; 

10. From Eq. ( p.l|) one can see that N{'K,bt,y) < 1. Therefore, saturation of the parton 
density means that 

N{x,bt,y) — > 1 at X — > 0; (2.2) 

11. Parameter k which was introduces in many papers ( or parameter W of Ref. [|l| ) is 
equal to 

«:(x,y) =Wi^,y) = N{^,bt = 0,yy, (2.3) 

12. The relation between parton density N{'x,bt,y) and gluon structure function xGa{x, ■\) is 
given by i i 

Ar(x,6, = 0,2/) = K{^,y) = X xGa{x,^), (2.4) 

which holds only for small N. We would like to stress that we do not need to know the gluon 
structure function for 1 since itself has a clear meaning as the scattering amplitude 

for a dipole. It should be also stressed that Eq. ( p.4| ) shows that A^(x, bt = 0, y) is not 
a parton ( colour dipole ) density p = xGa{x, ^)/7rR\ but it is the scattering amplitude 
which is equal to adipoie x P; 

13. Widely used function 0(A;^, y) (see Eq.(2.18) of Ref. for example ) is equal to momentum 
image of A^, namely, 

j d%t iV(x,6i,i/) = I d^k{l- e''^-'^ } (t){k\y) ; (2.5) 

14. as is used for as = asNc/n; 

15. Ra is the radius of a nucleus in the Gaussian parameterization of the nucleon density. In 
this parameterization the profile function in bt - representation looks as 

1 -4 

SA{bt) = — ^ e ^ . (2.6) 

R\{Gaussian) = ^R\{Wood — Saxon), where R^{Wood — Saxon) = r^A^ with tq = 
1.3 /m (see Ref. p for more details ); 
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16. 



All physical obsevables, related to nucleus, are denoted with subscript A while the nucleon 
observables will be marked by subscript A^. Both of these subscripts could be missed if the 
meaning of the physical quantity is obvious; 



17. As has been shown ^ the interaction between partons leads to a new scale Qcr(^? ^t) = 
1/x^j. for hdQCD. The line r = ln(l/x^) = \nQlj.{x,bt) is called critical line ( see Fig.l 
) and all physical quantities defined or calculated on the critical line are denoted with the 
subscript cr . For example, we denote a new scale, which has meaning of the average parton 
transverse momentum in hdQCD parton cascade, by Qlri^y h) while in other papers it has 
a different notation: Qlix,bt) = Ql{x)^^, = Q^x)^ = fi\x)'§^; 

18. All physical quantities defined to the right of the critical line ( for r > rcr = lnQ^^(x, 
will carry subscript — 

19. All physical quantities defined to the left of the critical line ( for r < = hiQl^{x,bt) ) 
will carry subscript 



2.2 The nonlinear equation for hdQCD evolution 

Yu. Kovchegov in his proof 0, that the GLR equation is able to describe the whole kinematic 
region including very low values of x, uses heavily two principle ideas suggested by A. Mueller 

m-- 

• The QCD interaction at high energy does not change the transverse separation between 
quark and antiquark ( the colour dipole size ), and, therefore, colour dipoles can be consid- 
ered as correct degrees of freedom at high energies which diagonalize the strong interaction 
matrix; 

• The process of interaction of a dipole with the target has two clear stages: 

1. Decay of the dipole into two dipoles, which is described by 

|^(xoi ^ xo2 + xi2)|' = ; (2.7) 

Xo2 Xi2 

2. Interaction of each dipole with the target with amplitude N{-x,bt,y). 
The equation is pictured in Fig. 2 and it has the following analytic form: 

N{^A,y) + ^ / rf^X2^^x (2.8) 

Vr Jp Xq2X-^2 



dN{yioi,bt,y) 
dy 



2 Cf as / Xn 



TT 
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y = ln(l/x) 

N>1 a3«l 

Quasi-classical 

Gluon Fields 




N= 1 



- N(r,b^,y)/d y dr = 
2agCj/7i N(r,bj,y) 
{ 1 +dN(r,b^,y)/dr}^ 

N(r,b^,y) = 
I^dr' N(r',b^,: 



n(r,bj,y) = N(r,bj,y) e' 

d n(r,bj,y)/d y dr = 
a^Cpg/Ti { 2 n(r,bj,y) 



Shadowing Corrections 



>r = Inl/xoi 



Figure 1: Different regions for the solution to the hdQCD evolution equation in the kinematic plot of 
DIS. The equation N{x,bt,y) = 1 gives the critical line. 



6 



2 A^(Xo2,bt - ^Xi2,?/) - A^(Xo2,bt - ^Xi2,l/) A^(Xi2,bt - ^Xo2,Z/) 

Assuming that Xi2 and X20 are much smaller than the values of the typical impact parameter ( 
bt ), one can see that Eq. ( p.8| ) is the GLR equation but in the space representation. It turns out 
that this representation is very convenient for searching solutions to Eq. ( |2.8|) . 




Figure 2: Pictorial representation of the nonlinear evolution equation. 

The first term in the r.h.s. of the equation gives the contribution of virtual corrections, which 
appear in the equation as a result of the normalization of the partonic wave function of the 
fast colour dipole ( see Ref. ). The second term describes the decay of the colour dipole 
with the size xqi into two dipoles wth sizes X02 and X12 and their interactions with the target 
in the impulse approximation ( notice factor 2 in Eq. ( [2. 81 ) ). The third term corresponds to 
simultaneous interaction of two produced colour dipoles with the target and describes the Glauber- 
type corrections for scattering of these dipoles. 
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Being differential equation, Eq. ( p.8|) needs an initial condition at ?/ = to be solved. In Ref. 
Qit was shown that the initial condition for N{'Xoi,bt,y) is the Mueller - Glauber formula ||15|| , 
namely, 

iV(xoi,6„i/ = 0) = 2(1- exp{--^^x^iAxoG^^^^^(a:o,4-)e"'*}l ! (2-9) 

y 4:I\cT^ri^ Xoi J 

where xoG^'^^^^(xq, -2-) is the gluon structure function obtained as a solution to the DGLAP 

"■01 

linear evolution equations M. 



2.3 The strategy of searching for solutions 

Eq. ( |2.8| ) is nonlinear integro-differential equation which is rather difficult to solve. However, this 
equation can be simplified and can be reduced to nonlinear but defferential equation in partial 
derivatives in two different kinematic regions (see Fig.l ): to the right of the critical line (r > rcr) 
and to the left of the critical line (r < r^r)- 



2.3.1 r > rcr 

Indeed, for r > rcr, the size of the two produced colour dipoles ( X02 and X12 ) are much larger 
than the size of the initial colour dipole (xqi ). Therefore, we can reduce the kernel of Eq. ( ^.8[ ) 
toi 

.i^x,^ x?../*^ (2.10) 

^02 ^12 "'^oi I ^02 ) 

Introducing a new function n(xoi,6t,|/) = iV(xoi, 64, ?/)/xqj^ and using the fact that the virtual 
corrections does not contribute in the double log approximation ( see Ref. |p for example ), one 
can obtain the equation: 

= ^ {2rUr,k,v) - e-'nUr,Ky)}; (2.11) 

ay dr tt 

where r = ln(l/xQj^). 

In section 3 we will discuss the solution to this equation which satisfies the initial condition of 



Eq. (P 



2.3.2 



r < Tc 



As have been discussed many times ( see for example Ref. |T3| and Ref. [|14| ), the new scale 
appears in hdQCD which has a simple physical meaning of the average transverse momentum 
of parton in the parton cascade ( Qcr{x,bt) ). In this region the size of the initial colour dipole 
is larger than the typical size which we expect in hdQCD parton cascade ( Xq-^ > 1/Qcr{x,bt) 
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). Therefore, the main contribution in Eq. (|2.8|) comes from the configuration when one of the 
produced colour dipole has a size much smaller than the size of the initial colour dipole Q We 
anticipate that the size of the smallest colour dipole will be of the order of 1/Qcr{x,bt). Such a 
configuration simplify the kernel in Eq. (|2.8|) which has a form 



c/'x2 3 TT / + TT / ^ ; (2.12) 

l/Qcr{x,bt) ^02-^12 Jl/Qcr{x,bt) Xq2 Jl/Qcr{x,bt) ^12 

The sum of two terms reflects the fact that two different produced colour dipoles can be small. In 
momentum representation Eq. ( p.l2| ) means that we sum log"(Q^j,(a;, bt)/Q^) which are the normal 



contributions for the DGLAP evolution coming from integration over transverse momentum from 
the small transverse momentum (Q^ ) to the large one ( Ql^{x,ht) ). 

Introducing a new function iV(r, 6j, y) = dr'N{r', bt, y), we reduce Eq. ( p.8|) to the differential 
equation: 



d^N^jr, bt,y) 2asCF r ,7 . , \ , m ( h ^ dN^r, bt,y) 

dydr = ^^{N^(rA,y) + N^{rA,y) }■ (2.13) 



One can see that in differential equation ( Eq. ( p.l3| ) ) we lost any dependence of the scale 



Qcr{x, bt) and it will come back to the problem only in matching a solution to Eq. ( ^.13] ) with a 



solution to Eq. (|2.11 



In section 4 we will find a solution to this equation which match the solution of Eq. ( p^.llD on the 
critical line. 

2.3.3 The general solution 

Let us summarize what solution we are looking for: 



1. We are going to solve the differential equation (see Eq. ( p.ll|) ) to the right of the critical 
line which satisfies the initial condition given by the Mueller-Glauber formula (see Eq. (|2T 



2. We are going to find the solution to the differential equation ( see Eq. (|2.13|) ) to the left of 
the critical line; 



3. We match the solution of Eq. ( 2.11|) with the solution of Eq. ( |2.13|) on the critical line; 



4. We can hope that the solution of these two differential equations will be close to the solution 
to Eq. ( p. 81 ) in the vicinity of the critical line (r ~ rcr)since we provide the matching of 
these two solutions on the critical line. 



-'^This idea has been suggested in Ref. in one of the first version of this paper but disappeared in the final 



version. 



9 



2.4 Experience of solving the GLR-type nonlinear equations 



During the past two decades we have learned the main properties of the solution to the GLR-like 
nonlinear differential equations. Our experience in solution of these equations based on three 
different approaches that have been developed. 

First, these equations were solved in semiclassical approximation, assuming N = where S* is a 
smooth function of r and y, namely dPS/dydr <^ {dS/dyY or (dS/dr)"^ (see Refs. [OlIU 
and references therein). In this approximation the GLR - type equation can be solved using 
characteristics method which leads to existence of a special ( critical ) line. All charecteristics in 
the region of low x cannot cross this line but can approach it. The nonlinear term provides that the 
parton density N{r, bt = 0, y) is constant on critical line. To the right of the critical line {r > r^r 
) the nonlinear corrections turns out to be small and the solution can be found as a solution 
of the linear DGLAP evolution equation but with the boundary conditions: = Const on the 
critical line. For our attempts to find a solution the important message from the semiclassical 
approach is the fact that we have a new scale {Qcr{x, h) ) and the properties of the solution looks 
differently for r > r^r = lnQ'^{x, bf) and for r < r^r- 



In Ref. |]rT|| a new method was suggested for solving Eq. ( |2.11| ). The generating function was 
constructed which linearize Eq. reducing Eq. ( |2.11| ) to a linear equation in partial derivatives 

but with one more variable. It is easy to find a general solution of this linear equation and all 
difficulties are concentrated in finding the solution which satisfies the boundary condition. We 
are going to use this method to find a solution to the right of the critical line. 

In the region of very low x, a new idea for searching a solution of the GLR-type nonlinear 
equation was suggested in Ref. fT^ . In this paper it was argued that the solution of Eq. (|2.11| ), 



which satisfies all physics restrictions , is a function of one variable: z = 4:asy — r. We will find 
this solution in an explicit way in section 4 and will show how to match this solution with the 
solution to the right of the critical line. 



2.5 Theory status of the equations 

In this subsection we recall the main assumptions that have been made to obtain the nonlinear 
evolution equationsfor hdQCD (see Eq. Eq. (|2J"T1) and Eq. (pl3D ). 



1. For derivation of all equations we used the leading log(l/x ) approximation of pQCD. Only 
in this approximation we can neglect the energy recoiled during emission of the colour dipole 
( the linear term in equations ) or during rescattring of this dipole ( the nonlinear term ). 
Formally speaking, we have to assume that parameter tts' log(l/x) > 1 while <C 1; 



2. To obtain the differential equations ( see Eq. ( [^.IID and Eq. ( [^.13| ) ) we have to assume 



existence of a more restricting parameter a^ry > 1, but a^r < 1 and asy < 1, which 
corresponds to the double log approximation of pQCD in the linear evolution equation; 
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3. The number of colours should be large ( A^^^ ^ 1). It has been shown 1 18 1 that for iV, 



the correlations between colour dipoles should be taken into account which break the simple 
physical picture of Fig. 2. In the double log approximation these correlations can be taken 
but they lead to rather complicated equations. On the other hand, in 
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mto account 

Ref. |jll| it was demonstrated that all these complications lead to corrections of the order 
of 1/iV^ and can be neglected even for Nc = 3; 



4. The master equation ( see Eq. ( p.8| ) ) sums so called "fan" diagrams ( see Fig. 3a ), which 
reflect the fact that the fast colour dipole decays into two colour dipoles. However, the 
recombination of two colour dipoles into one have been omitted (see Fig. 3b). For r > 
it has been proven that such diagrams give only small corrections |jl| but in the kinematic 
region to the left of the critical line such enhanced diagrams can be neglected only in the 
case of scattering off the heavy nucleus. It turns out ( see Refs. [|19| p ) that in the 
case of heavy nucleus the contributions of "fan" diagrams are proportionsal to 7A3 while 
enhanced diagrams of Fig. 3b are of the order of 7 without large factor A3. Therefore, 
our approach is valid only for heavy nuclei and can be considered only as a model for the 
interaction with a hadron; 



5. We assume that (i) there are no correlations between different nucleons in a nucleus and (ii) 
the average bt for colour dipole - nucleon interaction is much smaller than Ra- Both these 
assumptions are usual for treatment of nucleus scattering. 




Fig. 3-a Fig.3-b 

Figure 3: "Fan" ( Fig. 3a ) and enhanced ( Fig. 3b ) diagrams. 
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3 Solution to the right of the critical line 



In this kinematic region we use the generating function method, proposed in Ref. |TT|, which 
hnearize Eq. ( |2.11|) . Following Ref. |jTl[], we introduce a generating function 



9Air,y,bt,ri) 



oo 

E 

n=l 



(3.1) 



where n^-* = n^{r,y,bt) = N{r,ht,y) = {asu"^ /2Nc) p{r,bt,y) is proportional to the density 
of colour dipoles in the target, is proportional to the probability to have n - colour dipoles 
with the same size in the parton cascade. In the case, when we neglect correlations between 
colour dipoles in the parton cascade, = Cn{A)(^n^:^^ . Coefficients C„ (A) should be found 
from the initial condition of Eq. ( |2.9| ). 

For n^\r,bt,y) we can write the evolution equation of the Eq. (|2.11| ) - type, noticing that every 
colour dipole with size r at rapidity y and impact parameter bt can either propagate to rapidity 
y + dy changing its size from r to r + dr { see first two terms in Fig. 2 ) or decay into two colour 
dipoles ( see the last term in Fig. 2 and Fig. 3a). This observation leads to equation [|TT 



d'^n^_^\r,bt,y) - r 2 j, \ -r (n+i)/ 7 

— — = as{n n'i,>{r,bt,y) - ne n\J > {r,bt,y) } , 

oyor 



(3.2) 



where coefficient a^n^ corresponds to uj^yni^) where 7„ is the anomalous dimension of rSj^ and u 
is variable conjugated to y. In Refs. ||2CI|| the limit of small uj ( low x ) for 7„ has been calculated. 

Comparing Eq. (|3.1|) with Eq. ( |2.9|) one can obtain the following information on the generating 
function of Eq. (Ol) : 



(-1)" 



N{r, bt,y) 



2" n! ' 

2 QA^r, y, bt, r] = 



-r + r]A, 



n 



in) 
N,- 



n 



{r,y = 0{x = Xo), bt) = C„ (n;^^^^(r, y = 0(x = Xq), 

xoG'7v(xo,e'') S^ibt) ; 



DGLAP 



{r,y = 0{x = Xo),bt) 



2NcTTKi 



VA 



In { A R%SA{bt 



(3.3) 
(3.4) 
(3.5) 
(3.6) 
(3.7) 



Eq. (|3.5| ) and Eq. (|3.6| ) mean that we consider a; = xq to be so small that we can neglect correlations 
between produced colour dipoles. 

Eq. ( p.2|) can be rewritten as a linear equation for the generating function of Eq. (|3.1| ): 



dydr 



9A{r,y,bt,T], 



(3.8) 



as 



dr]' 



■9Air,y,bt,'n) - as'je 



■r—r) 



dgA{r,y,bt,r]) 
dr] 



9Air,y,bt,r]) 
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with 7 = 1. 

This equation is a hnear equation which can be solved just going to Melhn transform with respect 
to y and ^ = r + rj: 

gA{^,y,h,r^) = I ^^g^(^,co,h,p)e-y^^\ (3.9) 

In Eq. ( p. 91) the contours of integration over u and p he along the imaginary axis to the right of 
all singularities in u and p. 

Substituting Eq. ( p.9|) into Eq. (|3.8|) we obtain the following equation for the Mellin image 
gA{^,uj,bt,p). 

Tc = — \P - 1{P-^) e } gA{^,uj,bt,p) , 3.10 

which has an obvious solution: 

gA{^,u;,h,p) = ^?^(c.,p,6,)e^^^« + ^(P-^)(-^-^"'-^) . (3.11) 



Function gA{yj,p, ht) should be found from the initial condition of Eq. (|2.9| ) ( see also Eq. ( p.3|) - 
Eq. ( p.6| ) ). Our statement is that 

gA{.A,p) = n-fip)) ^e^^^-^^^-^^e^^^- (3.12) 

satisfies initial conditions, with f{p) = p + Substituting Eq. ( p.l2| ) in Eq. ( p.ll| ) and 

changing variables of integration from p to f ,we obtain 

9.K.V.,.) = /|^r(-/) + (3.13) 



where p{f) is determined by 



as 



f = P{f) + —pV) ■ (3.14) 



Indeed, aX y = { x = xq ) , \^\ <ti r, ^ = tja, r ^ 1, we have 



gAirM,y = yo,v = -r + VA) = I r(-/)e 



/ —nr I "5" I nrii / ijinn ^ 



E 



CO ( 1 An 



i:^(^x§,^..G--(..,-L)) (3.15) 
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Taking into account Eq. (|3.4| ), one can see that Eq. (|3.15| ) reproduces correctly the initial condition 
of Eq. (pD. 



It should be stressed that the contour of integration over / was taken along positive real axis in 
such a way that all positive and integer / are enveloped by it. The integral over / is well defined 
due to rapid decrease of r(— /) ( see Ref. [|lT| ). To study asymptotic behaviour of the parton 
densities at low x we have to move the contour toward the imaginary axis. The asymptotic stems 
from the rightmost singularities in / as well as from possible saddle points in this integral. The 
first dangerous region is a vicinity of / ^ 1 where we expect a saddle point. To investigate a 
situation near f= 1, we introduce a new variable f = 1 + t, hoping that t will be very small. 
At t — > our solution of Eq. ( p. 131) has a form 



d^^^i.,r,tM ^VA+v ^ (3.16) 



where f] = In (]^). 



In Appendix A we calculate the integral of Eq. (|3.16|) using the saddle point approximation and 
the answer, that we obtain in this appendix ( see Eq. ( |A.5| ) ) can be written in the form: 



,2 ^,_^2/l 



9A,Ay, e = VA, r, h) = J^7 e «i e--^-^+^ , (3.17) 



with ojcr is given by Eq. (|A.3| ) and Iq by Eq. (|A.4| ) 



One can see that Eq. ( p.l7| ) leads to QAiVi i = f]Ai^i W) which is constant on the critical line given 
by equation: 

^ = UcrVcr — T + r + r]A_ + f] — lu ^ — - = Const{as) (3.18) 

where Const{as) can be a function of as but it does not depend on neither y nor r or rjA- We 
anticipate that the value of to will be small on the critical line and, because of this, the last term 



in Eq. (|3.18|) has been taken into account in the phase 



The solution of Eq. (|3.18| ) gives 



Vcr = -^(r - 2{r,A + f,) + 2). (3.19) 

Aas 



This line gives the minimal value for namely, \E' = 1. 
The value of gA{y,C, = Va, r, ht) is equal to 



gA,cr{y^^ = VA,r,bt) = (3.20) 

TT 
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Solving Eq. ( 3.19|) with respect to r we obtain the typical momentum ( colour dipole size ) on 



the critical line 



At 6i = Qlri^] bt = 0) oc As, as can be see easily from Eq. (|3.21| ). 

Comparing Eq. (|3.2CI|) with Eq. ( p^) one sees that xGa{x, ^) on the critical line is equal to 

xGaA^A) = ^^-^-^^/^- iV(x^6, = 0,y,,) = ^R\Ql{x) oc aI . (3.22) 

It is interesting to notice that such A-dependence can be understood directly from the Mueller - 
Glauber formula of Eq. (|2.9| ). Indeed, the shadowing corrections become to be essential when 

2 

-x^ AxGix, l/x^) ^ 1 . (3.23) 



2 



Recalling that anomalous dimension in the vicinity of the critical line is equal to 7(0;) = = 1) 
one can find from Eq. ( p.23|) that x^^. = l/Q^r °^ A~^ , which is the same A-dependence that 



we obtained from more advanced calculations ( see Eq. (|3.22|) ). 



4 Solution to the left of the critical line 



In this section we are going to solve Eq. ( |2.13| ) using the idea suggested in Ref. |0, namely, 
= N^^{z) where 

z = Aasy - r - l3{A-ht) . (4.1) 

Function (3{A] bt) can be chosen in a such way that the equation for the critical line in terms of 
new variable z looks as z = 0. One can see from Eq. ( p.l9|) that j3{bt) is equal to 



4 E?2 9/,2 

Pih) = -2^ - 2 1n^ + ^ + 2- (4.2) 
Assuming that N^{r, bt, y) is a function of z and bt only, we can rewrite Eq. (|2.13|) in the form: 

Here we used the explicit equation for variable z, namely, z = AasV — r — (3{bt). 
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4.1 A solution to Eq. ( |^) 

First, we find a solution to Eq. (|4.3| ) introducing function ({z) in the following way: 

N^{z,bt) = dz' (^l - e-^(^'**)) . (4.4) 



For function C,{z) Eq. ( |4.3| ) can be reduced to 

dazM) 



dz 

Changing variable in the integral from z' to Q' we have 



^ f dz' [l - e^«(^''^*)) . (4.5) 



dz 4 Jcoibt) ^ 



ill (^--n- 



Eq. ( [4.6|) can be easily solved and the solution is 

, = V2 r , . (4.7) 

J(o{bt) + (e-C - 1) 

The numerical solution to Eq. ( [4.7] ) is given in Fig.4. This solution depends on the initial condi- 
tions on the critical line . There are two of them, namely. 



N^{z = 0) = N^{z = 0) ; (4.^ 

InJV^I _ dlnN^ l _ 1 //i ( 

dz 1^=0 — dz 1^=0 ~ 2 ■ 



We can assume that ( ^ 1 at large In this case we can neglect e in the integrand and we 
obtain an explicit analytic solution: 



or in terms N{z, bt) we have 



N^{z,bt) = 1 - e~^v^+^^' (4.11) 

where zo{bt) should be found from matching of this solution with the solution of the previous 
section on the critical line. 



4.2 Matching with the solution to the right of the critical line 



We cannot use the simple solution of Eq. ( [4.10| ) for matching with the solution of Eq. ( p.l7| ) 



since C is expected to be rather small in the vicinity of the critical line and, therefore, we cannot 
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neglect e in Eq. ( [4.7|) . However, for small ( <^ 1 we easily obtain a simple solution to Eq. (|4.71 ), 
namely, 

^ - (4.12) 



In- 



z . 



Co{bt) 2 

Eq. ( [4.12|) leads to Eq. ( [4.9|) of matching condition. This equation allows us to find Coih), namely, 

1 - e-^o(''*) = Ncr,^ . (4.13) 



Eq. ( ^.12| ) and Eq. ( [4.13| ) provide matching as well as determine the parameters of the asymptotic 
behaviour for given by Eq. (|4.11|) . Using Eq. ( p.20| ) one can find from Eq. ( [4.13| ) that 



Uh) = 2 



eas 



TT 



(4.14) 



N(z) 




2 . 5 



10 



12.5 



20 
Z 



Figure 4: Dipole number density N{z) as a function of the critical line parameter z. Solid line is the 
numerical solution with N'(0) = 0.01, while dashed line corresponds to N'{0) = 0.4. The bold 
dashed line shows the asymptotic solution of Eq. with N'{0) = 0.4. 



Fig. 4 illustrates the behaviour of N{z) versus z. One can see that asymptotic is reached only at 
very large z ~ 7 — 10 or for very low x ( x ~ 10^^ — 10^'' ). At HERA we have z ^ 4 — 5 and 
we are far away from the asymptotic solution. However, we can penetrate the region of large z 
using a nuclear target, as one can see in Eq. ( |3.19| ). Indeed, for heavy nuclei we can easily have 
2riA ~ 3.5 — 4 . It gives a possibility to have even at HERA kinematic region sufficiently large 
values of z ~ 7.5 — 9. 
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4.3 Stability of the solution 



Arguments of the last section show, that the solution given by Eq. ( ^4.7| ) has a big chance to be 
the solution to our equation. To prove this statement we need only to show that the solution 
of Eq. ( [4.7| ) is stable with respect to small perturbations of the initial conditions on the critical 
line. In other words, let us consider a small variation of the initial conditions \6(cr{y,'r,bt)\ < S 
where < 5 <C 1. The solution is stable if for every 6 we can found a small e such that e ^ 
at 6 0. We can find a linear equation for function 6({y,r,bt) substituting ({y,r,bt) — > 



Cs + ^CiVi^ibt) where C,s is given by Eq. (p?7|). Assuming that 6(cr{y,r,bt) is small we obtain a 
linear equation 

= as Say,r,b,) e-<^^^K (4.15) 

ay or 

In variables z and r this equation has a form 



d'^6C{z,r, bt 



1 = 1 SC{z,r,bt) e-^=(") , (4.16) 



d z dr 4 

which can be solved going to Mellin transform with respect to variable r. Indeed, for Mellin image 
we have 

-^5C(^,.,6,)e-^=(^). (4.17) 



dz Av 
The general solution of Eq. ( [4.17| ) looks as follows 



ic(., 6,) = /°"°° ^ m^v - i . (4.18) 

J a—ioo Z TT t 

The integral over z' in Eq. ( [4.18| ) does not depend on z at large z. It means that the solution 
given by this equation is actually the function of r only. It is obvious that such a solution cannot 
be tolerated by the initial conditions on the critical line. 

4.4 Low X ( high energy ) asymptotic 

As we have discussed the total dipole cross section is equal to 

C^dipole = (4.19) 

//I ^2 „ OO /* ^ 

d\2N{r,y,bt) = 2n { ° db^^ + db^^ } = 2n j db^ = 27ibl{y,r) 
Jo Jbl Jo 

since N^{r, y, bt) — > 1 at low x ( see Eq. ( ^.11| ) ) and N^(t, y, bt) gives a small contribution to 
the integral in Eq. ([4.19|) . 

The value for bl{y,r) we can evaluate recalling that solution of Eq. ( |4.11| ) is written for 
Aasy -r + 2fi + 2\n^ - ^ > 0. Therefore, 



z 



bl{y,r) = ^ (Aasy-r + 2\n^] , (4.20) 
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which gives 



71 R\ I iasy-r + 2\n 



For the gluon structure function Eq. ( [4.21| ) leads to 



xG{x,Q'^ 



Q'R\ 4 as ln(l/x) - InQ^ + 2 In 



N 



R\ 



(4.21) 



(4.22) 



Eq. ( [4.19D can be rewritten as an integral over z, namely 



//•OS roo 
d%2N{r,y,bt) = 2n { dhj + / dhj N^ ; 
Jo Jbl 

+ 2aseR\ 



Aasy-r + 2\n-Ji- 



dz 1 



(4.23) 
(4.24) 



It is interesting to notice that cr^ipoie for nucleus target manifests itself a remarkable scaling in the 
kinematic region to the left of the critical line. Indeed, from Eq. ( [4.24| ) one can conclude that 



^dipole 



a 



N 

dipole 



Xb 



7lR% 



(4.25) 



Fig. 5 shows the behaviour of the ratio 



dipole 



at a; = 10 ^ with y = ln(0.01/x) for different 



nuclei. The full line describes the nucleon target and other lines correspond to nuclei with A= 
40 (Ca), 150, 296 (gold) going from bottom to the top. One can see that the cross sections are 
quite larger than the geometrical estimates. 

It should be stressed that Eq. ([4.20|) is an artifact of the oversimplified Gaussian parameterization 
for the nucleus profile function ( see Eq. ( |2.(j| ). In more realistic approach to the nucleus profile 
function ( for example the Wood-Saxon one ) instead of expl—b'f / R\) in Eq. (|2.6| ) we have a 
function which is equal to 1 for bt < Ra and falls down as exp(— 6^//;,), where h does not depend 
on A, at ht > Ra- Such function gives bo{y, r) 

h 



AD2 

bo{y,r) = Ra + - \Aasy-r + 2\n^ 



(4.26) 



instead of Eq. (|CT|)Pi 

Therefore, we expect the following asymptotic behaviour for adipoie and xG{x, Q'^) 



'^dipole 
2\ 



f h ^ AR 
271 \R + -{Aasy - r +2ln^ 



N 



} 



xG{x, Q 



\R + --{Aasy - r+2\n^] 



h 
2 



2 \ 2 
i 

I 

Ri 



(4.27) 
(4.28) 



We are very grateful to Al Mueller for pointing us a model feature of Eq. (4.2C). 
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Figure 5: The ratio of the dipole cross section to the geometrical estimates ^'^^'^ for low x = 10 ^ 

( y = ln(0.01/x) ) as a function of for different nuclei. All other curves correspond to 
A=40, 150, 300 going from the bottom to the top; as = 0.25, = IGeV"^ . 



5 Summary 

In this paper we found the solution to the evolution equation for high parton density QCD ( see 
Eq. ( |2.8| ) ). This solution is given by two equations: Eq. ( |3.13| ) to the right of the critical line ( 
see Fig.l) and Eq. ( |4.4|) to the left of the critical line described by Eq. (|3.19| ). 



This solution gives the colour dipole density (p = xGa{x, ■\)/itR\ ) which is small to the right 
of the critical line , reaches the value of the order of unity on the critical line and increases up to 
the value p ~ l/a^ at very small values of x. 

Such a behaviour of p reveals itself in the following properties of the gluon structure function: 



1. To the right of the critical line xG is given by the DGLAP evolution equations; 



2. On the critical line xG{x,Qlj.{x)) = Ql^.{x)R\ (see Eq. (|3.22|) ) where ,Ql^{x) is defined 
by Eq. (CT) ; 



3. To the left of the critical line xG{x, Q"^) = Q"^ R\{Aasy - r + 2 In 



As far as A-dependence is concerned, in the kinematic region to the right of the critical line we 
have xGa oc AxGn, while on the critical line xGa{x,QI^{x)) oc As and only to the left of the 
critical line we have xGa oc Q^R\ oc As. 
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Our solution reproduces the saturation of the gluon density but due to sufficiently mild dependence 
on the impact parameter the saturation leads to the dipole-target total cross section proportional 
to ln(l/x) in the region of the extremely low x (see Eq. ( [4 .211) for the exact answer). 

Both the exact solution near the saturation limit (see Eq. ( |4.11| ) ) and the proportionality to 
xG (X Q^R\ ( see Eq. ([4.22|) ) are the manifestation of the fact that nonlinear equation leads to 
the parton distribution which peaks at the size 1/Q^^,[x). In other words, in the region of 

low X the parton distribution has a definite mean transverse momentum, namely, < >= Q'^j.{x). 
Therefore, this solution supports the more intuitive approach that has been discussed before for 



the parton distribution at low x [113] 



We want also to draw your attention to the fact that in the region of small x our solution leads to 
a new scaling between DIS with nucleon and DIS with nucleus given by Eq. ( |4.25| ). This equation 

2 

allows us to calculate how nucleus DIS approaches A3 dependence in the region of low x. 

We hope, that our solution will stimulate a more detailed study of the properties of the high 
density parton system and, in particular, it will lead to more quantitive development of ideas 



suggested in Ref|]T3 



A Appendix 

In this appendix we calculate the phase ip{uj, r, t, bt) in Eq. ( ]3.16| ) at t ^ 1. Using Eq. p.l4| ), we 
can find p for 05 <^ 1 and t ^ 



p = 1 + t - — (l + 2t) , (A.l) 



which leads to 



il:{uj,r,t,bt) + ti7]A + fi) = e "^"^{p - I) + pr + ujy + t{'qA + r]) 

UJ 

^ Ljy + r - —r - tr I 6-"^-^ - 1 + + -{va + v)] ■ (A.2) 

u \uj{—r) ijj r j 

Integration over t gives a 5 - function which defines the value of a; = iOcr- 



where ^ = 1 + 



We can make the integration over u = Ucr + A , considering ^ -C 1. 

ilj{uj,r,t,bt) 

I a or 2aQtr ^ 



asr ( asr 2astr\ 
-r + V uJcrV + A I y 
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-r H h uJcrV + A- 



UJr. 



V 



y 

V 

to 



(A.4) 



Integration over A gives 5 ^"g( ^) (^^ _ which yields 



9Aiy,^ = VA,r, ht) 



Co-; 



2Ti2as{r)tQ<, An R\ 



e A e ' 



(A.5) 
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